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abstract 

Radar  backs  catter  from  a  chaff  cloud  is  examined  for  the  case  in 
which  the  radar  approaches  the  cloud  from  above  and  moves  into  and 
through  the  cloud.  The  chaff  cloud  is  modelled  as  a  collection  of  randomly 
distributed  and  randomly  oriented  dipolea  having  a  mean  volume  density 
that  is  a  function  of  altitude  but  invariant  in  a  horizontal  plane.  The 
statistics  of  the  chaff  echo  power  when  die  number  of  dipoles  in  the 
resolution  volume  is  large  are  reviewed.  The  impulse  response,  which 
approximates  the  short* -pulse  radar  response,  and  the  radar  return  signal 
for  a  specific  range  resolution  cell  as  a  function  of  the  position  of  the  radar 
relative  to  the  cloud  ar  i  derived  for  several  chaff  cloud  configurations  and 
antenna  patterns.  Pov/er  spectral  density  functions  for  range-gated  chaff- 
return  signals  and  correlation  functions  for  the  signals  following  square- 
law  detection  are  also  derived. 
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1.  INTRODUCTION 


Radar  systems  are  frequently  required  to  operate  in  the  presence  of  echoes  from  chaff, 
ground  and  sea  surfaces,  clouds,  and  the  various  forms  of  precipitation.  To  operate  effec¬ 
tively  in  such  environments,  the  radar  must  be  designed  to  minimize  the  effects  of  these 
clutter  signals  which  may  lead  to  spurious  detections,  that  is,  false  alarms,  or  to  desensitiza¬ 
tion  of  the  radar  so  that  the  actual  target  is  shielded.  A  prerequisite  for  radar  design  optimi¬ 
zation  is  a  model  for  the  clutter  environment  that  is  sufficiently  general  to  cover  the  range  of 
conditions  likely  to  be  encountered  while  retaining  the  simplicity  essential  to  its  utility  as  a 
design  instrument. 

This  report  examines  a  model  for  a  chaff  cloud  consisting  of  a  collection  of  randomly 
oriented  and  randomly  distributed  dipoles.  Specifically  considered  is  the  response  of  an  air¬ 
borne  radar  transmitting  a  short  pulse  and  moving  rapidly  relative  to  the  chaff  cloud.  The 
impulse  response  function  for  chaff  backs catter  is  derived,  and  the  return  from  a  fixed  range 
is  examined  for  a  radar  moving  through  the  chaff.  The  results  obtained  are  valid  for  narrow 
transmitted  pulse  widths  and  may  be  extended  to  cases  involving  wide  pulses  by  convolving 
the  transmitted  pulse  envelope  with  the  impulse  response  function. 

The  first-  and  second-order  probability  density  functions  for  chaff  echo  powers  are  re¬ 
viewed,  and  the  correlation  and  spectral  density  functions  associated  with  this  chaff  model  are 
derived. 


2.  ENCOUNTER  REFERENCE  FRAME 


For  the  following  discussion,  it  is  assumed  that  the  chaff  cloud  is  confined  to  a  specified 
region  of  space.  In  all  cases  the  chaff  cloud  will  be  assumed  bounded  above  by  a  plane  surface 
parallel  to  the  tangent  plane  to  the  mean  surface  of  file  earth.  A  Cartesian  reference  frame  is 
constructed  with  the  radar  at  the  origin  as  shown  in  figure  1. 

It  is  further  assumed  that  the  radar  velocity  is  constant  throughout  the  encounter  which 
implies  that  the  flight  angle  y  between  the  z-axis  and  the  radar  velocity  vector  is  also  constant. 

The  radar  receives  returns  from  chaff  elements  located  in  that  portion  of  the  spherical 
shell  lying  within  the  cloud  and  with  boundaries  at  radial  distances  d  and  d  +  Ad  from  the  ra¬ 
dar.  The  distance,  h,  from  the  radar  to  the  upper  edge  of  the  cloud  is  po  itive  when  the  radar 
is  above  the  cloud  top  and  negative  when  below.  For  a  radar  speed  V,  this  distance  varies 
during  the  encounter  as 


h  =  d  -  V  t  cos  y,  (1) 

with  the  time  origin  taken  as  the  time  when  the  radar  range  gate  first  penetrates  the  cloud  to 
a  depth  equal  to  the  range  gate  width.  Ad. 

3.  CHAFF  CLOUD  MODEL 

The  chaff  cloud  is  assumed  to  consist  of  a  collection  of  resonant  dipoles,  randomly  orien¬ 
ted  and  randomly  distributed  within  a  volume,  large  relative  to  the  radar  resolution  cell.  It 
is  further  assumed  that  the  dipole  locations  are  statistically  independent,  so  that  the  number 

Preceiing  page  blank 
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Figure  1.  Encounter  reference  frome.  (The  distance  h  from  me  rodor  to  the 
top  of  the  chaff  cloud  is  token  as  positive  when  the  rador  is  above  the  cloud 
top  and  negative  when  below). 


of  dipoles  found  in  any  subvolume  is  a  Poisson  distributed  random  v.-\riable.  It  should  be  noted 
that  such  a  distribution  does  not  account  for  the  effect  of  chaff  chimj  ing,  or  "birdnesting." 

Since we  assume  that  the  basic  scatterer  of  the  chaff  cloud  is  a  dipole,  it  is  appropriate 
to  lode  first  at  the  scattering  behaviov  of  a  single  dipole.  The  radar  backs catter  from  a  di¬ 
pole  is  usually  maximum  when  the  dipois  is  viewed  at  broadside  with  the  electric  vector  of  the 
incident  wave  parallel  to  the  axis  of  the  dipole.  The  magnitude  of  this  return  as  a  function  of 
the  length  i  of  the  dipole  is  shown  in  figure  2.  At  l/\  _  0.45,  the  curve  shows  a  definite  peak; 
Because  of  this  resonant  cross  section  enhancement,  half-wave  dipoles  give  the  'maximum  re¬ 
turn  per  unit  length  of  the  wire.  Dipole  thickness  has  a  relatively  small  effect  od  the  return; 
the  general  tendency  being  that  the  amount  of  resonant  peaking  is  diminished  for  thicker  dipoles. 

As  a  dipole  is  rotated  in  the  plane  formed  by  the  direction  of  incidence  and  the  electric 
vector  of  the  incident  wave,  the  amplitude  of  the  backs cattered  field  changes.  In  figures  3 
and  4,  the  radar  cross  sections  are  sbowD  as  functions  of  rotation  for  a  half-wave  resonant 
dipole  (i/\  -  0.452)  and  for  a  dipole  whose  length  is  2.46  \  ,  respectively.1  It  may  be  ob¬ 
served  that  in  the  first  case  the  cross  section  exhibits  a  single  lobe;  but  for  the  longer  dipole. 


S.  ond  Liepo,  V.,  ’Measured  Bode  Scattering  Cross  Section  of  Thin  Wires,'  The  Uni-crsity  of  Michigan 
Radiation  Laboratory,  Report  No.  8077-4-T,  1967. 
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°/a  4 


Figure  2.  Backscattering  cross  section  of  a  wire  par  nil  el  to  electric  vector;  a/\- 6.27  x  10~^.  (a  s  radius 
o*  wire,  't  -  length  of  wire,  \  rwovelength.) 

multiple  lobes  appear.  For  certain  lengths,  the  amplitude  of  the  sidelobes  may  be  greater 
than  that  of  the  broadside  lobe.  It  may  also  be  noted  that  in  the  case  of  the  half-wa/e  dipole 
the  |sins  4>'\  curve  provides  a  close  empirical  approximation  to  the  experimental  curve 
(fig.  3). 

If  the  dipole  is  rotated  in  a  plane  perpendicular  to  the  direction  of  incidence,  a  similar 
variation  in  radar  cross  section  occurs.  In  this  case  the  cross  section  has  been  found  to  vary 
as  sin4  S',  where  the  angle  S'  is  defined  in  figure  5.  Hence,  we  have  the  half-wave  dipole 
radar  cross  section  as  a  function  of  the  polarization  oi  the  incident  field  (horizontal  as 
depicted  in  picture  fig.  5)  and  the  dipole  orientation  along  the  two  principal  planet  of 
rotation. 


Accordingly,  we  can  write 


'if!*?  =0.8*  !sin*  <*>'!. 

o  '=t  n 


-  (6*.  <*>') 

^  *  '=  IT. 


0.89  sin4  0‘. 


^deBettoncourt,  J.  T.,  "Bistntie  Cross  Sections  of  Cylindricol  Wires,*  Pickard  &  Bums,  Inc.,  Scientific  Report  No.  1, 
P&B  Pub.  No.  735A,  Wolthom,  Mass.,  1961. 


X  '  - 


j  CENTER  LINE 
I  OF  WIRE 


WIRE 
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Figure  5.  Geometry  for  computation  of  the  cross  section 
of  a  dipole. 

Assuming  that  a (6 ',d>')  =  0.89  \7  ct,  (<*>')  a2  ($,')  where  c?1  (&')=  sin4  S'  and  cr7  (<£')  = 

|  sin5  4-‘  I  which  seems  to  be  a  valid  approximation  in  the  case  of  the  half-wave  dipole, 
we  may  combine  equations  (2)  and  (3)  to  obtain  an  empirical  angular  expression  for  a 
half-wave  resonant  dipole  as 

^  ^  ^  =  0.89  sin4  S'  jsin5<p'|.  (4) 

\2 

Since  half-wave  dipoles  are  most  efficient  in  that  they  provide  the  maximum  return  per 
unit  weight,  we  shall  be  concerned  hereafter  with  such  dipole  lengths  only. 

In  a  ctnff  cloud  consisting  of  a  large  number  of  dipoles,  it  is  impossible  to  specify  the 
orientation  of  ~.h  of  the  individual  dipoles.  We  are  therefore  confronted  with  the  problem  of 
determining  an  average  dipole  cross  section  based  on  die  probability  density  of  dipole  orien¬ 
tations.  We  consider  two  possible  distributions:  (a)  uniformly  distributed  over  4n  steradians, 
and  (b)  confined  to  a  horizontal  plane  but  uniformly  distributed  over  Zn  radians. 

In  case  (a),  the  probability  density  of  dipole  orientations  becomes 


p  (&',  <t>')  -  _L  sin  S' 

4  77 


0  S.  4>'  <  2  77 


0  s.  S'  <  t; 


The  mean  cross  section  per  dipole  is  independent  of  incident  polarization  and  direction  of  in¬ 
cidence  and  is  obtained  simply 
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where  V  =  V  (8',  4>‘  ,8  ,<P)iB  a  random  function  of  the  direction  of  incidence  (8,  4>)  and  the 
dipole  orientation  angles  (8',4>‘)  with  respect  to  a  spherical  reference  frame  centered  on  the 
dipole  (figs.  6  and  7).  The  rotation  angle  4*  is  measured  in  a  horizontal  plane. 


Considering  first  a  horizontally  polarized  incident  field,  we  obtain  from  figure  6  the 
following  relations 


Z  HZ 


Figur*  6.  Gtometry  fcr  horizontal  potori ration  calculations. 


=  !  tan*y 


Sin 


sin2  0  +  tan2  ¥{ 


5/2 


Hence,  we  have 


?h  (*>  .  0.89  f 

X2  ‘  2  77  J» 


JTT 


(1  -  cos2  ¥  cos2  £)2 


tan2  ¥ 


sin2  &  +  tan2  ¥ 


5/2 


d¥. 


<1< 


which  nan  be  evaluated  numerically  (see  fig.  8  for  plot)  or  by  direct  integration  at  &  -  0,  tt/2 
to  obtain 


oH  (0)  =  0.334  X2. 


(i: 


5,  (77/2)  =  0.302  X2. 


<i: 


Figure  8.  Mean  backscolter  cross  section  as  a  function  of  elevation  [&)  of 
incident  wave-horizontal  polarization. 


When  the  incident  field  is  vertically  polarized,  the  random  variable  H*  may  be  seen  to  co¬ 
incide  with  the  dipole  elevation  angle  S'.  Thus,  we  obtain 


X2 


0.89  jcos*  sj 
_ 


s  i  n4  <p  dT 


(13) 


=  0.334  |coss  S\ 


As  to  the  actual  distribution  of  dipole  orientations,  the  available  evidence  appears  to  sup¬ 
port  the  horizontal  distribution  at  least  for  X-band  dipoles.3  However,  enoatraospberic  dis¬ 
persal,  the  selective  weighting  of  dipoles  to  assume  other  attitudes,  the  recent  seeding  of 
dipoles,  and  the  buffeting  effect  of  turbulent  air  currents  may  give  rise  to  chaff  clouds  better 


'Palermo,  C.  J.  and  Boyer,  L.  H.,  "Bistatic  Scattering  Cross  Section  of  Chaff  Dipoles  with  Application  to  Com- 
muni  cations,*  Proceedings  IEEE,  Aug.  1965- 


described  by  the  uniform  distribution  of  dipole  orientations.  In  any  event  the  mean  radar 
cross  sections  per  dipole  for  both  the  horizontal  and  uniform  orientations  are  employed  in 
subsequent  calculations  involving  impulse  and  range-gated  radar  responses  to  various  chaff 
clouds. 

To  consider  the  radar  return  from  each  of  the  dipoles  contained  within  an  arbitrary 
volume,  we  must  first  define  a  chaff  dipole  volume  density.  As  noted  previously  we 
assume  that  the  number  of  dipoles  falling  within  any  given  region  v  is  a  Poisson  dis¬ 
tributed  random  variable,  Nv .  We  shall  parameterize  Nv  with  respect  to  space  coor¬ 
dinates  -  that  is,  define  a  random  process  in  space  by  identifying  Nv  with  each  point 
P(Xq,  y0,  z0)  as  follows.  Let  vp  be  the  region  defined  with  respect  to  a  Cartesian  frame 


V 

a 


f.  .  Ax,,  Ax 

(x.  y.  z):  x0-t*x£x0  +  _. 


Ay 

2 


*y  *y0 


:  s  2  s  z. 


and  the  volume  encompassed  by  this  region  is 


Av=AxAyAz 


Now,  let  N(vp)  be  the  number  of  dipoles  contained  within  vp .  We  assume  that  the  mean  value 
of  N(vp),  ElN(vp)},  and  the  limit  cl  the  ratio  of  the  mean  value  of  N(v  )  to  the  volume  Av  as 
the  volume  shrinks  to  the  point  P  exist  and  are  continuous  functions  of  the  space  coordinates. 
We  define 

m  (vp>  ;E(N  (vp)h  (14) 


and 


S(P) 


A 


i  im 


Ax-0 


m  (v  > 
v  p' 

A  v 


Ay  -0 
A  z  -  0 


(l*r') 


where  S  (P)  is  the  desired  chaff  dipole  volume  density.  It  should  be  noted  that  our  definition  of 
chaff  volume  density  avoids  the  usual  difficulties  of  attempting  to  take  tne  limit  of  a  discrete 
process  (see,  for  example,  Moore*).  Although  N(Vp )  is  discrete  so  that^lim  N(v  )/Av  cannot 

exist,  this  is  not  the  case  for  lim  E{N(v  )}/Av).  p  0 

Av-0  p 


The  choice  for  the  chaff  dipole  density  function  presents  a  problem.  A  realistic  model 
for  this  function  is  not  available;  it  depends  on  such  parameters  as  the  type  of  chaff  used, 
method  of  dispersing,  and  prevailing  winds.  Our  analysis  is  confined  to  two  basic  models:  an 
exponentially  increasing,  semi-infinite  distribution;  and  a  uni  modal  distribution  which  initially 
increases  to  a  maximum  value  and  then  gradually  decreases.  Constant-density  surfaces  in 
both  cases  are  horizontal  planes  with  the  zero  density  plane  at  a  distance  h  from  the 
radar. 


^Moore,  R.  K.  and  Williams,  C.  S.,  "Rador  Return  at  Near-Verticol  Incidence,*  Proceedings  IRE,  Vol.  45, 
1957. 
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The  first  model,  hereafter  referred  to  as  the  exponential  model,  is  defined  by 


S  (h')  HTT~tl  -  e-a<h-h'>]  h  £  h' 

=  0  h  >  h' 

where 

h  =  distance  from  the  radar  to  the  chaff  cloud  top, 

5  (h‘)  =  the  chaff  dipole  volume  density  at  the  distance  h*  from  the  radar, 
j  =  maximum  volume  density, 

o 

a  =  a  parameter  related  to  the  magnitude  of  the  density  gradient,  that  is 


lv  S  (h'V=.h  I  =  a  s 


It  should  be  observed  that  in  the  limit  as  ,  the  normalized  exponential  density  function, 
S  (h  ‘)/&a,  becomes  the  Heaviside  unit  step  function.  Hence,  this  model  also  serves  to  describe 
the  uniform-density  chaff  cloud  with  an  infinite  density  gradient  at  h '  =  h. 

The  other  model,  hereafter  referred  to  as  the  Rayleigh  model  because  of  its  similarity  to 
the  Rayleigh  probability  density  function,  is  defined  as  follows: 


s  (h')  =  n  <x(h' -  h)  exp  |--|(h' -  h)2J  hsh'  (17) 

=  0  h  >  h' 

where 

n  is  the  total  number  of  chaff  dipoles  contained  in  the  vertical  column  with  a  unit 
err.-.'-  ~?d?onal  area. 

a  is  a  parameter  related  to  the  density  gradient  and  maximum  volume  density  5a  by 
the  relationships 


IV  5  (h')h.=h!  =  n  a 


(18) 


Sa  =  ?(b')!h'.hM/^=n|/?.  (19) 

An  important  limiting  case  for  this  model  is  obtained  again  by  allowing  a-*.  The  result 
in  this  case  is  the  Dirac  delta  function  with  intensity  n.  Hence,  this  model  in  the  limit  de¬ 
scribes  a  surface  layer  of  dipoles. 


Although  both  models  as  detailed  above  are  assumed  to  have  infinite  extent  along  the 
horizontal,  such  an  assumption  is  not  critical  to  the  analysis  to  follow.  In  reality,  no  chaff 
cloud  possesses  either  the  uniformity  or  the  infinite  extent  of  the  models  postulated.  It  is  a 
simple  matter  to  eliminate  the  infinite  extent  criterion  by  merely  nutting  off  the  cloud  at  some 
distance  from  the  radar  greater  than  the  maximum  range  cell  incorporated  in  the  radar  under 
consideration.  Whenever  the  chaff  cloud  has  dimensions  that  are  large  relative  to  that  raoge, 
the  cloud’s  effect  on  the  radar  will  be  essentially  as  modelled.  Regarding  the  criterion  of 
horizontal  uniformity,  here  again  we  are  concerned  from  a  practical  standpoint  with  uniformity 
over  that  region  seen  by  the  radar.  If  the  extent  of  the  cloud  is  large  relative  to  the  observed 
region,  then  the  nonuniformity  due  to  tailing  off  along  the  periphery  can  be  discounted. 

From  a  computational  standpoint,  the  use  of  the  Rayleigh  model  is  particularly  convenient. 
The  parame!  jr  n  alone  determines  the  number  of  dipoles  contained  in  the  chaff  cloud.  Thus,  a 
could  be  alloyed  to  change  as  the  cloud  evolves,  but  so  long  as  n  is  constant  the  total  number 
of  dipoles  remains  constant.  No  such  conservation  rule  is  applicjible  to  the  exponential  model 
with  the  possible  exception  of  the  infinite  depth  cloud,  in  which  case  an  unlimited  number  of 
dipoles  is  involved. 

As  noted,  a  relates  to  the  density  gradient  and  the  dispersion  of  the  chaff  cloud.  It  there¬ 
fore  constitutes  an  inverse  age  description  for  the  cloud.  Thus,  a  will  be  large  immediately 
after  dispersal  (high  maximum  density  gradient)  and  will  decrease  as  the  cloud  settles  (lower 
maximum  density  gradient:).  Since  a-  is  also  rolated  to  the  location  of  peak  density  in  the  Ray¬ 
leigh  model,  ffie  value  at  any  given  time  r  after  dispersal  may  be  determined  from  the  approxi¬ 
mation 


a  =  (vch  T>'2-  (2°) 

where  Vq,  is  the  mean  fell  rate  for  the  chaff  dipoles. 

The  feet  that  a  may  be  a  function  of  time  introduces  nonstationarity  and  consequently  non- 
ergodicity.  We  shall,  however,  be  concerned  with  the  properties  of  the  chaff  cloud  during 
measurement  intervals  that  are  very  short  when  compared  with  the  reciprocal  of  the  time 
derivative  of -Cn  a.  On  this  basis,  we  shall  assume  that  the  process  is  quasi-stationary  in  the 
sense  that  its  statistical  properties  are  essentially  invariant  during  the  measurement  internals. 

This  paper  disregards  two  effects  common  to  volume-distributed  scatterers:  multiple 
scattering  and  shadowing.  Since  both  effects  become  important  in  dense  distributions,  we  shall 
establish  a  simple  criterion  for  the  degree  of  rarity  of  scatterers  such  that  both  effects  are 
negligible.  Multiple  scattering  arises  as  a  consequence  of  the  fact  that  the  reradiated  field 
from  each  dipole  is  a  linear  superposition  of  reradiatioo  of  the  primary  incident  field  as  well 
as  reradiation  of  the  fields  scattered  by  surrounding  dipoles.  Shadowing  occurs  when  the  pri¬ 
mary  field  is  attenuated  in  passing  through  intervening  layers  of  scatterers  before  illuminat¬ 
ing  file  deeper  layers  of  dipoles.  Both  of  these  effects  become  prominent  if  the  totality  of  the 
scattered  fields  is  of  the  same  order  of  magnitude  as  the  primary  field.  If  such  is  not  the 
case  that  is,  if  the  primary  field  predominates  throughout  the  scattering  region  then  quite  ob¬ 
viously  multiple  scattering  can  be  disregarded.  Further,  if  little  energy  is  removed  from  the 
primary  field  and  reradiated  as  a  scattered  field  (and  little  energy  is  dissipated  thermally 
within  the  scatterers),  then  the  attenuation  due  to  the  presence  of  the  scatterers  of  the  primary 
field  as  it  traverses  the  scattering  region  may  be  neglected. 

Since  the  fraction  of  the  energy  removed  from  the  primary  field  and  backs cattered  toward 
the  radar  is  a  good  indication  for  randomly  oriented  scatterers  of  the  relative  magnitude  of  the 
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primary  ami  scattered  fields,  we  can  obtain  an  upper  limit  on  this  ratio  by  considering  the 
sum  of  scatter er  effective  cross  sections.  That  is,  we  can  formulate  our  criterion  by  requir¬ 
ing  that  the  sum  of  the  solid  angles  subtended  by  the  scatter  er  cross  sections  be  small  over 
any  angular  region  considered.  Hence,  we  obtain  the  requirement 


*k  4.  4  77 


o+A*  r*.  - 

4, 


dd 


(•9q+&9 

I  |  sin  8  dS  d^> 

Jg.  Jd, 


<  <  1, 


(21) 


where 

B,  is  the  maximum  distance  range  cell; 

Ro  is  the  minimum  distance  from  the  radar  to  the  scatterers;  angular  intervals  (£0,  $0  + 
A  9)  and  (<£0,  4>0  +  M>)  are  arbitrary. 

By  considering  angular  intervals  sufficiently  small  and  applying  the  mean  value  theorem,  we 
obtain  the  requirement 


■  90'  dr  <  <  1 


(22) 


If  S  (r,  d,  4)  is  monotonic  increasing  and  equation  (22)  holds  along  the  direction  of  steep¬ 
est  ascent  —  that  is,  the  direc  .on  determined  by  V  5  -  then  the  inequality  holds  everywhere. 
Therefore,  if  multiple  scattering  and  shadowing  may  be  disregarded  along  the  gradient,  they 
may  be  disregarded  everywhere. 

For  a  nonmonotonic  density  such  as  the  Rayleigh  cloud,  this  criterion  can  be  used  to  ex¬ 
amine  the  relevance  of  multiple  scattering  and  shadowing  in  the  vertical  direction  for  Rk  =  =. 
Inserting  equation  (17)  into  equation  (22)  and  integrating,  we  obtain  the  results 


o  n/4  77  <  <  1  hiRp, 

_  -o/jfjn-h)1 

one  ^  /4  77  <  <  1  ^  >  h  i  0, 

_  -a/ 2  af 

o-  n  e  _/4  77  <  <  l  h  <  0. 


(23a) 

(23b) 

(23c) 


Taking  o  =  0.33  K2  and  using  criterion  (23a),  we  have 


n  <  <  37.74A2- 


(24) 
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Considering  X-band  dipoles  (\  -  .05  m),  we  find  the  result 


n  <  <  15,000/m2 


as  a  sufficient  condition  for  disregarding  multiple  scattering  and  shadowing  effects  along  the 
vertical.  Recognizing  that  n  is  the  number  of  dipoles  over  a  square  meter,  we  see  that  such 
effects  usually  are  totally  negligible. 


4.  MEAN  1MPUISE  RESPONSE 

4.1  Signal  Return  from  a  Single  Scatterer 

We  consider  first  a  radar  transmitting  a  pulse  P{  (t)  toward  a  target  with  radar  cross 
section  o-  and  located  at  a  distance  d  from  the  radar.  The  backs cattered  received  power  is 
then  given  by 

I  i  a  G.  G,  K--  o 


/  2  d\  a  “t  *•"  2  d 

Pr  (t)  =  Pt  t  -  ±1) - L_! - t  i  ±1 

\  c  /  (4  it)3  d4  c 


t<l£ 


wher'.  Gt  and  Gr  are  the  gains  in  the  target  direction  of  the  transmitting  and  receiving  anten¬ 
nas,  respectively;  d  is  the  distance  from  the  radar  to  the  scatterer;  and  c  is  the  speed  of  light. 
For  simplicity  we  shall  assume  that  Gt  =  Gr  =  G;  hence,  Gt  Gr  =  G2. 

Strictly  speaking,  equation  (26)  is  correct  only  for  continuous  wave  transmission,  since 
the  width  of  the  frequency  spectrum  is  nonzero  in  every  other  ca?,  and  therefore  the  wave¬ 
length  \  is  not  a  constant.  In  addition  the  radar  cross  section  and  antenna  gains  are  also 
functions  of  the  wavelength.  Nonetheless,  whenever  the  spectrum  of  the  transmitted  pulse  is 
narrow  band,  this  equation  provides  a  valid  approximation  to  the  backs  cattered  power  if  we 
insert  for  k  the  wavelength  corresponding  to  the  center  frequency. 

If  the  target  considered  is  randomly  oriented,  the  received  power  is  also  random.  We 
find  its  mean  value  by  taking  the  expectation  of  equation  (26)  to  obtain 


p,  (t) = pt  (t  -  —)  aGly 

'  C  •  (4  77)3  d4 


4.2  Signal  Return  from  Volume  Distributed  Scattererg 

We  assume,  as  stated  previously,  that  the  individual  scatterer  locations  are  statistically 
independent.  Further,  we  shall  assume  that  the  volume  to  be  considered  has  a  depth  in  range 

that  is  large  relative  to  the  wavelength,  and  that  the  scatterer  volume  density  is  nearly 
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constant  over  a  distance  large  relative  to  the  wavelength.  Under  these  constraints,  the  co¬ 
herent  return  from  the  volume  will  be  entirely  negligible  when  compared  with  the  incoherent 
return.5  We  shall  therefore  concentrate  on  the  incoherent  return  only. 

Employing  the  dipole  volume  density  function  defined  previously,  we  now  define  a  chaff 
volume  cross  section,  o-° , 


(r.  e ,  <t>)  -  or ^  S  (r,  0,  #), 


(28) 


where  the  subscript  v  refers  to  the  common  transmit  and  receive  polarization  state.  (Al¬ 
though  no  major  modification  is  required,  we  shall  not  treat  cross-polarized  transmission- 
reception  states  in  this  report.) 

We  shall  consider  the  signal  return  from  a  volume  AV,  containing,  by  hypothesis,  N(AV) 
individual  dipoles.  As  before,  N(AV)  is  a  Poisson  distributed  random  process.  The  incoherent 
radar  backs catter  from  such  a  volume  is  obtained  by  summing  the  contribution  of  each  scat- 
terer  contained  in  the  incremental  volume.  If  the  incremental  volume  subtends  a  small  angu¬ 
lar  sector  so  that  the  antenna  gain  may  be  considered  constant  across  the  sector,  the  back- 
scattered  power  Is 


APr  (i)  = 


G2  (<?.  <t>)  K2 
(4  n)3 


N(  Av) 


(29) 


where  d  is  the  distance  from  the  radar  to  the  fth  scatterer.  We  assume  that  the  dimensions 

l 

of  the  incremental  volume  are  small  relative  to  the  mean  distance  to  the  radar;  that  is 
dj  >  >  3/A  V,  then 


APr  (t)  =- 


G^2p«  {-—)  ^  ^ 
(i17)S  /— ~ 


(30) 


where  explicit  dependence  on  the  angles  0  and#  has  been  suppressed.  Assuming  the  a,  to  be 
distributed  with  mean  5^  and  with  m(AV)  as  the  mean  value  of  N(AV),  we  obtain  for  the  ex¬ 
pected  value  of  APr  (t) 

_  G2  \2  Pt  (t  - 

APf(t)= - 1 - L'  m(AV)5  (31) 

(4  v)3 


A.  and  Gold  twin,  H.,  'Coherent  ond  Incoherent  Scattering  from  Assemblies  of  Scatterers,'  in  Propaga¬ 
tion  of  Short  Radio  Waves,  Kerr,  D.  E.,  ed.,  McGraw-Hill,  New  York,  1951. 
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Dividing  both  sides  by  AV,  taking  the  limit  as  AV  shrinks  to  the  point  (d, &,4>)  and  using  equa¬ 
tion  (15),  we  obtain 


—  G2  A2  p 
dp(t)  _ 

d  v 


■•(-¥) 


(4  tt)3  d4 


S  (d.  9,  4> )  cr 


Substituting  equation  (28),  we  have 


dpr  (t) 


G2  A2  P, 


2pt  (t-^) 

- - - -  °°v  (d.  <t>) 

(4  v )3  d4 


The  mean  received  power  scattered  from  an  arbitrary  volume  V  is  then  given  by 


Pr  (t)  =  _M_  f  Pt  ft  d  v 

(4  17)3  J  *  \  C  /  <j4 


where  and  G  are  functions  of  location  and  aspect,  respectively.  Equation  (34)  is  a  general 
expression  for  the  mean  backs  carter  from  chaff  but  is  not  especially  useful  in  its  present 
form.  To  carry  the  analysis  further,  we  shall  reduce  this  expression  to  a  more  tractable 
form,  but  let  us  first  clarify  the  use  of  equation  (34)  by  considering  several  specific  examples. 


;  from  a  Semi-Infinite  Cloud 


Consider  the  geometry  shown  in  figure  9.  An  isotropic  antanna  of  gain  G  is  located  at  a 
height  h  above  the  chaff  cloud  which  has  uniform  >.nsity  and  unbounded  depth.  A  spherical 
coordinate  system  with  its  origin  at  the  radar  is  used.  In  such  a  case,  equation  (34)  becomes 


g2<a2  r  r  r2w  pt  (*  “4-) 

?r  (t)  = -  I  j  |  - - - Lt2  sin 

(4*)3  i  Je  Jo 


6  d$>  d<?  d  r. 


where  8  =  cos-1  (-h/r)  and  we  have  assumed  a  constant  so  that  it  has  been  removed  from 
the  integral.  Integrating  equation  (35)  in  4>  and  6 ,  we  have 


P,  (t) 


G2  o°  A2 

"n 

2  (4  7T)2 


fp,  (t-— )-  --] 
A  '  c  'U*  rJJ 


Introducing  a  change  of  variables  2r/c  =  t  ,  and  using  the  fact  that  Pt  (t  -  2r/c)  =  0  for 
1 5  2r/c,  we  obtain 

G2a°\2  I*  ri  it,  ti 


G2a°A2  i*  ri 

(t)= - —  P,(t-T)  1-Hll  dr. 

(4  7t)2  c  Jj h/c  j2  C  t* 


ANTENNA 


Figure  9-  Cross  section  of  o  semi-infinite  chuff  cloud. 


Equation  (37)  can  now  be  recognized  as  a  superposition  integral  with  Pt  (t)  as  the  driving  func¬ 
tion  and 

(4  77)2  c  (t2  c  t3J 


=  0. 


t  > 


2  h 


t  £ 


2  h 


(38) 


aa  tbs  response  to  an  impulse  function.  We  should  note  that  equation  (38)  is  not  formally  cor¬ 
rect,  since  by  postulating  an  impulsive  transmitted  waveform  we  have  violated  our  premise 
on  the  extent  of  the  illuminated  region.  However,  the  impulse  response  function  does  provide 
a  valuable  tool,  since  a  solution  may  be  obtained  by  convolving  the  impulse  response  with  any 
transmitted  waveform  that  has  sufficient  spatial  extent  so  that  this  premise  is  not  violated. 
Further,  convolution  with  a  transmitted  waveform  with  a  spatial  extent  of  several  wavelengths 
will  not,  in  general,  provide  a  result  significantly  different  from  the  impulse  response  func¬ 
tion.  Hence,  we  may  consider  the  impulse  response  function  as  an  approximate  solution  to  the 
chaff  backs catter  problem  whenever  the  transmitted  pulse  has  a  duration  very  short  as  com¬ 
pared  with  the  impulse  response  but  a  spatial  extent  of  at  least  several  wavelengths. 

If  we  introduce  another  variable  d  =  ct/2,  which  may  be  interpreted  as  the  distance  from 
which  the  signal  is  arriving  at  time  t,  we  have 


second  part  is  a  function  of  time  and  specifies  the  shape  of  the  returned  pulse. 
4.2.2  Finite  Thickness  Chaff  Layers 


The  return  for  a  finite  thickness  uniform  density  chaff  cloud  follows  directly  from  equa¬ 
tion  (39).  If  the  upper  edge  of  the  chaff  layer  is  at  h  a  and  the  lower  edge  at  h2,  then 


P;(t)  =  0, 


(40a) 


P,  (t)  =- 


G2  \2 
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There  is  the  regular  range  delay  until  the  pulse  hits  the  cloud.  Then  before  the  pulse 
reaches  the  trailing  edge  at  h2,  the  return  is  the  same  as  that  from  a  semi-infinite  cloud;  mid 
after  the  leading  edge  of  the  pulse  has  passed  through  the  chaff  layer,  the  return  is  given  by 
equation  (40c). 


4.2.3  Varying  Density  Chaff 


We  can  find  the  differential  power  return  from  a  chaff  layer  of  infinitisimal  thickness  by 
allowing  (h2  -  hj),  equation  (40c),  to  shrink  to  dh  giving 


dP‘  (t)  = 


G2  \2 


(4  77)2 


_ 2 

c  t 


i.  dh. 
2  d 


(41) 


The  volume  cross  section  cr°  ,  which  up  to  now  has  been  treated  as  a  constant,  can  new 


be  treated  as  a  variable  when  we  integrate  both  sides  of  equation  (41).  The  returned  power 
from  a  varying  density  chaff  cloud  is  then 
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(42) 


Psr  (t) 


g2\2 
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Having  examined  some  simple  examples  and  introduced  the  concept  of  the  chaff  impulse 
response,  let  us  return  to  the  more  general  form  of  equation  (34).  We  note  that  for  those 
dipole  orientation  distributions  considered  in  section  3,  3^  has  the  form  k2  Cn  |cos  9\  k(’’>  . 
(In  the  horizontally  polarized,  horizontally  distributed  case  we  may  use  the  approximation 
v  =  0.32 X2.)  We  shall  consider  hypothetical  antenna  gain  functions  of  the  form  G(0)  = 

G0  j  cos  9  i  .  Equation  (34)  becomes 


Go  c» 

P  (t)  =---  1— 
(4  77)3 


r  r  r  p‘  (*-¥) , 

Jh.  4  Jo  r* 


cos  6\  2^+k  sin  9  S  (r,  6)  d(pd9  d  r.  (43) 


where 


5  (r,  0)  =  Sffl  (1  -  exp  (-  a  r  (cos  8  -  cos  9))),  8  i  9 


=  0, 


e  <  9 


for  the  exponential  chaff  model,  and 

f  (r.  8)  -  n  a  r  (cos  6  -  cos  9)  exp  (cos  9  ~  cos  9)j.  9i9 

=  0.  e  <  9 


for  the  Tlayleigh  chaff  model,  with 


Q  =  < 


0  |h/r|  >  1 

cos-1  (-h/r)  jh/r|  <1,  and 


hn  =  max  (0,  h). 

Integrating  with  respect  to  4>  and  using  an  impulse  driving  function,  we  have 


(44) 


(45) 


(^Jj,  |c°*  6\litk  sin9  S  fSJ.  ej  d  B. 
Returning  to  our  previous  notation,  that  is,  h  *  =  -ct/2  cos  9 ,  we  obtain 


(46) 


p;  «>  I"'"  s  (h.  v,  „h-. 

4  (4  77)2  \c  t/  Jd 


(47) 


where  da  =  max  (h,  -ct/2).  Substituting  d  =  ct/2  and  defining  R(X  ,tj  )  =  cGjj  \«c  /4(4 a),  we 
obtain  ’ 
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. . . 


where 


J*  d  - 

h.k+2-e  h  (hi  h<)  dh  _  (SKr,d  )k  (  ’  h'**2*  S  (h,  h')  dh' 

Jo 


sgn  dm  = 


1  for  d  >0 


-1  for  d  <  0 


Upon  substituting  for  the  exponentially  varying  chaff  density,  we  obtain  an  integral  that  can 
be  reduced  by  successive  integration  by  parts  with  the  following  result: 


f  h*»  (1  -  e~tt<h'~h>)  dh'  =  -1--  +  s  !-a-  + 
Jn  s  +  1  a‘+» 


S  +  1  a«41 


d* 

Z_i  i  ! 


A  similar  application  of  successive  integration  by  parts  results  in  an  iterative  solution 
when  the  Rayleigh  chaff  density  is  substituted.  If  we  define 


=  f  h'*  (h'  - 
Jo 


h)  e‘a/2  <h  -h>  dh'. 


we  obtain  the  iterative  formula 


where 


»  x 

I  (S)  =  L,  +  s  hs_1  *  +  — I  (s  -  i  -  2). 


L  =  e~ah  n  -  e~a/2  (d*h)3  =  I  (0). 


=  -  d*  e~a/2  <d*h> 


s  =  1,  2, 


<*>=  ^  e~«/2  (h*-h)3  dh' 


We  shall  exemplify  the  use  of  these  formulas  by  considering  four  hypothetical  antenna 
patterns:  an  isotropic  pattern,  two  teardrop  lobe  structures,  and  a  hollow  cone  pattern  which 
may  be  represented,  respectively,  as 


Gj  =  Gf0  isotropic  pattern; 
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G|  =  G20  cos4  (0)  0  S  0  £  77,  double-lobe; 

G2  =  G^0  cos4  &  2-  £0  <77  single-lobe; 

0  Z9  <- 
2 

G4  =  4  Gj0  (1  -  cos2  6)  cos2  0  ~<  8  £  77  hoi  1  -w  cone, 

=  0  0  £0  <  — 

2 


(53) 


(54) 


(55) 


where  G10,  GJ0,  G30,  and  G40  are  the  peak  power  gains  of  the  assumed  antennas. 

We  obtain  the  chaff  impulse  response  for  a  radar  system  employing  these  antenna  when 
both  exponential  and  Rayleigh  cloud  models  are  considered  and  the  chaff  dipoles  are  uniformly 
oriented  for  both  polarization  states,  or  the  dipoles  are  horizontally  distributed  and  the  trans¬ 
mission-reception  polarization  is  horizontal.  We  also  consider  the  combination  of  an  isotropic 
antenna,  horizontally  distributed  chaff  orientations  and  a  vertically  polarized  wave.  The  spe¬ 
cific  impulse  response  functions  obtained  in  each  case  are  enumerated  in  appendix  A.  The 
results  are  also  plotted  as  figures  10  through  17. 

It  should  be  observed  that  the  impulse  response  functions  for  chaff  return  differ  signifi¬ 
cantly  from  those  obtained  from  ground  return.6  While  the  impulsive  response  from  surface- 
distributed  scatterers  obtains  its  maximum  value  initially,  the  impulsive  response  from  chaff 
rises  to  its  peak  value  only  after  the  incident  wave  has  penetrated  some  non-zero  depth  into 
the  scatter  volume.  The  amount  of  retardation  of  peak  return  relative  to  initial  return  depends 
upon  the  distance  from  the  radar  to  the  chaff  cloud  and  the  specific  model  considered.  In 
general,  the  impulse  response  for  chaff  peaks  most  rapidly  when  the  radar  antenna  i»  narrow- 
beam  and  the  magnitude  of  the  chaff  cloud  density  gradient  is  large.  It  may  further  _  t  ob¬ 
served  that  when  the  chaff  dipoles  assume  a  horizontal  orientation  and  the  polarization  is 
vertical,  the  net  effect  on  the  impulse  response  is  equivalent  to  the  effect  obtained  by  using  a 
narrower  beam  antenna;  that  is,  backs catter  from  the  vertical  direction  is  accentuated  where¬ 
as  returns  from  lateral  directions  are  strongly  attenuated.  Thus,  a  horizontal  orientation  of 
dipoles  with  vertical  polarization  increases  the  rapidity  with  which  the  impulse  response  peaks, 
with  the  consequence  that  the  chaff  cloud  response  more  closely  resembles  the  surface  scat¬ 
tering  of  rough  terrain. 

Of  perhaps  even  greater  importance  than  the  chaff  impulse  response  is  the  chaff  range¬ 
gated  response  that  is,  the  return  to  the  radar  from  a  specific  range  resolution  cell  extending 
from  d  to  d  +  Ad  as  the  radar  altitude  h  is  allowed  to  vary. 

The  mean  power  return  from  chaff  at  time,  t,  is  given  by  the  convolution  of  the  trans¬ 
mitted  pulse,  Pt,  with  the  mean  impulse  response 


%di«>n,  A.  R.,  Moore,  R.  K.,  and  Warner,  B.  0.,  "Radar  Return  or  Near-Vertical  Incidence,  Summary  Report,'  Univ. 
of  New  Mexico  Engineering  Experiment  Station,  Tech.  Rep.  EE-24,  Sept.  1959. 
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Figure  12.  Impulse  response  functions  at  h  =  100  m 
for  a  sinfcosd  antenna  pattern  and  uniformly  distrib¬ 
uted  dipole  orientations  (exponential  cloud  model) 


Figure  13.  impulse  response  functions  at  h=  100  m 
for  an  isotropic  ontrnno  pattern,  vertically  polarised 
wave,  ond  horizontally  distributed  dipole  orientations 
(exponential  cloud  model) 
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Figure  16.  Impulse  response  (unctions  at  h  =  100  m 
for  a  sinOcosd  antenna  pattern  and  uniformly  distrib¬ 
uted  dipole  orientations  (Rayleigh  cloud  model) 


Figure  17.  Impulse  response  functions  at  h  =■  100  m 
for  an  isotropic  antenna  pattern,  vertically  polarized 
wave,  and  horizontally  distributed  dipole  orientations 
(Rayleigh  cloud  model) 


pr  ft) 


J  Pt  (l  -  V)  Pr^>dr. 


(56) 


If  the  transmitted  wave  amplitude  is  constant  with  Pt  (t)  =  P0  ,  and  we  are  concerned  only 
with  the  return  contributed  by  dipoles  lying  at  ranges  between  d  and  d  +  Ad  from  the  radar,  we 
obtain  for  d  2  h 


A  Pr  (d)  = 


1 


d+Ad 


?0P, 


(r)  d  r. 


(57) 


which  becomes,  for  Ad,  sufficiently  short 


AP.  (d)  =  A  d  PQ  P*  (d). 


(58) 


When  rewritten  as  the  ratio 


APf  (d)/PQ  =  AdpS  (d). 


(59) 


we  obtain  the  chaff  backscatter  loop  attenuation  for  the  range  cell  (d,  d  +A  d). 

Considering  the  hypothetical  antenna  patterns  and  dipole  orientations  previously  intro¬ 
duced,  we  have  plotted  representative  chaff  backscatter  loop  attenuation  functions  versus 
radar  altitude  as  figures  18  through  27. 


5.  PROBABILITY  DENSITY  FUNCTIONS  FOR  CHAFF  BACKSCATTER 

Since,  by  hypothesis,  the  chaff  dipoles  are  randomly  and  independently  distributed  and 
randomly  oriented,  the  statistics  of  the  scattered  field  are  identical  with  those  for  the  sum  of 
sine  waves  with  random  identically  distributed  amplitudes  and  phases.  It  is  well  Imown  that 
such  a  process  is  asymptotically  Gaussian  when  the  number  of  constituent  waves  of  the  sum  is 
large.7  The  envelope  of  a  Guassian  process  will  be  Rayleigh  distributed  when  it  is  defined  for 
narrow-band  processes  in  the  manner  of  Rice'  and  when  the  broader  definition  of  envelope  as 
the  absolute  value  of  the  complex-valued  function  obtained  by  adding  to  the  given  function  i 
times  its  Hilbert  transform  is  used.  Thus,  we  should  expect  that  the  envelope  of  the  chaff  re¬ 
turn  process  will  be  Rayleigh  distributed  and  that  the  amplitude  of  the  squared  and  low-pass 
filtered  return  will  be  exponentially  distributed  for  a  sufficiently  large  number  of  chaff  dipoles 
contributing  to  the  return.8  Hence,  the  radar  cross  section  for  chaff  which  is  proportional  to 
tho  backs  catter  ed  newer  should  have  a  probability  distribution  which  is  asymptotically  expo¬ 
nential  as  the  number  of  contributing  dipoles  becomes  large. 


^Rice,  S.,  'Mathematical  Analysis  of  Random  Noise,'  Bell  Systems  Technical  Journal,  vol.  23,  1944  and  vol.  24, 
1945. 

®Borison,  S.  L.,  'Statistics  of  the  Rodar  Goss  Section  of  o  Volume  of  Chaff,'  MIT,  Lincoln  Laboratory,  Group  Rep 
1965-10,  1965. 


31 


t  H  /  > 

s*l  / 

§  *51  / 

o  r  / 

2  -  I 


ICO. CO  33.33  *33.33  -iGC.CO  *144.47  -233.33  -300.0 

h€!GKT  or  RADAR  ABOVE  CHAFF  CLOUD 

Figure  18.  Range-gated  response  at  d  =  100  m  for  an 
isotropic  antenna  pattern  and  uniformly  distributed 
dipole  orientations  (exponential  cloud  model) 
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Figure  20.  Ranae-goted  response  ot  d  =  100  m  for  a 
single-lobe  cos *0  ontenno  pottern  and  uniformly  dis¬ 
tributed  dipole  orientations  (exponential  cloud  model) 
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Figure  19.  Range-gated  response  at  d  =  100  m  for  a 
cos20  antenna  pattern  and  uniformly  distributed  dipole 
orientations  (exponential  cloud  model) 
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Figure  21.  Range-gated  response  at  d  •-  100  m  for  a 
sindcosd  ontenno  pattern  and  uniformly  distributed 
dipole  orientations  (exponential  cloud  model) 
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Figure  22.  Range-gated  response  at  d  =»  100  m  for  an 
isotropic  antenna  pattern,  vertically  polarized  wave, 
and  horizontally  distributed  dipole  orientations  (expo¬ 
nential  cloud  model) 
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Figure  23.  Range-gated  response  atd  =  100  m  for  an 
isotropic  antenna  pattern  and  uniformly  distributed 
dipole  orientations  (Rayleigh  cloud  model) 
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Figure  24.  Range-gated  response  at  d  =  100  m  for  a 
cos20  antenna  pattern  o.id  uniformly  distributed  dipole 
orientations  (Rayleigh  cloud  model) 
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Figure  25.  Range-gated  response  at  d  =  100  m  for  a 
single  lobe  cos'O  antenna  pattern  and  uniformly  dis¬ 
tributed  dipole  orientations  (Rayleigh  cloud  model) 
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Figure  2S.  Range-gated  response  at  d  =  100  m  for  a 
sindcosO  antenna  pattern  and  uniformly  distributed 
dipole  orientations  (Rayleigh  cloud  model) 
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Figure  27.  Range-gated  response  at  d  =  100  m  for  an 
isotropic  antenna  pattern,  vertically  polarized  wave, 
and  horizontally  distributed  dipole  orientations 
(Rayleigh  cloud  model) 


The  fact  that  the  instantaneous  voltage  variate  of  the  return  signal  is  Gaussian  in  the 
limit  is  a  consequence  of  the  central-limit  theorem  of  probability  theory  and  is  not  dependent 
upon  the  distribution  of  the  component  amplitudes  of  the  individual  dipole  returns-  The  rate 
b*.  which  this  limiting  form  is  approached,  however,  does  depend  on  this  distribution. 

Gouda mit  and  Schiff  have  examined  the  finite  number  of  scatterers  case  and  have  found  the 
Gaussian  distribution  to  be  a  good  approximation  when  ranges  to  the  contributing  dipoles  are 
all  equal,  even  when  the  number  of  dipoles  equals  4  or  5.  When  the  number  is  10  or  more, 
they  found  the  differences  are  completely  negligible.5  Borrison  has  also  examined  the  use  of 
the  exponential  law  for  describing  the  chaff  radar  cross  section  and  concluded  that  the  ap¬ 
proximation  is  good  when  the  average  number  of  dipoles  in  the  radar  resolution  volume  is 
greater  than  two  or  three.8 

Since  the  exponential  probability  law  is  a  single  parameter  distribution,  it  is  completely 
determined  by  its  mean  value.  Thus,  we  may  represent  the  first  order  distribution  for  chaff 
radar  m  oss  section  as 

P  (O  =  _L  e  0j>  /C\  (60) 


^Siegert,  A.  ond  Goldstein,  H„  "Coherent  ond  Incoherent  Scattering  from  Assemblies  of  Scatterers,*  in  Propaga¬ 
tion  of  Short  Radio  Waves,  Kerr,  D.  E.,  ed.,  McGraw-Hill,  New  York,  1951. 

®Borison,  S.  L.,  'Statistics  of  the  Rodor  Cross  Section  of  o  Volume  of  Choff,’  MIT,  Lincoln  Loborotory,  Group  Rep 
1965-10*  1965.’ 
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For  such  a  distribution,  we  note  that  the  variance  equals  the  square  of  the  mean,  and  that  in¬ 
cluding  more  dipoles  in  the  radar  resolution  cell  by  increasing  the  resolution  cell  size  does 
not  narrow  the  relative  distribution  width,  since  the  mean  to  rms  ratio  is  constant.  It  should 
be  observed  that  the  most  probable  value  of  backs cattered  power  is  zero.  The  return  power 
will  be  less  than  half  the  average  value  29  percent  of  the  time  and  greater  than  twice  the  aver¬ 
age  14  percent  of  the  time. 

As  the  radar  moves  relative  to  the  chaff  cloud,  the  path  lengths  to  the  individual  dipoles 
will  vary.  The  doppler  shifts  thus  introduced  will,  in  general,  be  different  for  each  dipole 
scatterer;  hence,  frequency  dispersion,  rather  than  mere  frequency  shift,  will  result  in  a  re¬ 
turn  signal  that  is  noise-like  and  narrow-band  when  the  transmitted  signal  frequency  is  greater 
than  the  maximum  doppler  shift.  Successive  samples  obtained  by  a  pulsed  radar  will  generally 
be  correlated,  the  degree  of  correlation  for  a  fixed  sampling  period  depending  on  the  width  of 
the  doppler  dispersion  spectrum. 

The  second  order  probability  density  for  sample  powers  Pj  and  with  sampling  interval 
r  becomes 


P2  (Pj.  P2,  r)  = 


^-(p, 'P2)/pao -p)  ^  ! 2  vp  Pj  P2^ 
P^d-P)  °\pa  (i-py 


(61) 


where  P.  is  the  average  returned  power,  and  the  correlation  coefficient  p{t)  is  defined  by 


P,  P,  -  P2 

p(r)=.  1  2  * 


(62) 


and  I0  is  the  zeroth  order  modified  Bessel  function  of  the  first  kind.  For  r  sufficiently  large 
so  that 


Pi  pi-prp2  =  pJ- 


(63) 


we  observe  that  p  (r)  -0;  and,  as  expected,  we  have 


P2  (Pj,  P2;  r)-p  (Pj)  p(P2).  (64) 

That  is,  when  the  sampling  interval  is  sufficiently  long  so  that  the  correlation  coefficient  is 
zero,  the  sampled  return  powers  are  not  only  uscorrelated  but  also  independent. 


6.  POWER  SPECTRAL  DENSITY 


Of  paramount  importance  in  determining  whether  the  return  powers  from  successive  re¬ 
turn  pulses  are  uncorrelated  is  the  width  of  the  power  spectral  density.  The  correlation  of 
detected  pulse  amplitude  from  pulse- to- pulse  affects  the  required  integration  time  to  achieve 
a  given  detection  probability;  for  high  detection  probabilities,  the  greater  the  correlation,  the 
longer  the  required  integration  time. 
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The  successive  range- gated  pulse  returns  from  chaff  may  be  considered  analogous  to  the 
sequence  of  sampled  values  obtained  by  periodically  sampling  a  random  process.  As  we  have 
previously  noted,  each  of  the  dipoles  may  be  thought  of  as  a  simple  harmonic  generator,  trans¬ 
mitting  with  random  phase  and  amplitude  and  with  a  frequency  determined  by  the  radar  trans¬ 
mission  frequency  and  the  doppler  shift  induced  by  its  motion  relative  to  the  radar.  The 
ensemble  of  returns  as  seen  at  the  radar  comprises  a  narrow-band,  noise-like  process,  the 
decorrelation  time  of  which  is  inversely  proportional  to  its  bandwidth.  If  this  decorrelation 
timo  is  short  compared  with  the  radar  pulse  repetition  period,  successive  pulses  will  provide 
independent  samples. 

Although  the  chaff  cloud  drift  velocity  and  dipole  rotations  are  major  contributors  to  the 
doppler  shift  for  a  stationary  radar,9  these  factors  are  not  significant  when  compared  with 
the  frequency  shift  induced  by  the  rapid  movement  of  the  airborne  radar  altimeter.  There¬ 
fore,  in  calculating  the  power  spectral  density,  we  shall  consider  only  the  radar  movement. 

In  calculating  the  chaff  return  power  spectral  density  we  note  that  the  doppler  frequency 
shift  is  the  same  for  signal  returns  from  dipoles  whose  position  vectors'  intersections  with 
the  radar  velocity  vector  form  equal  angles  (see  fig.  1).  The  amount  of  power  returned  at  any 
doppler-shifted  frequency  depends  therefore  on  the  antenna  gain  in  the  specified  direction,  and 
the  chaff  radar  cross  section  in  that  direction  and  at  the  range-gate  distance  from  the  radar. 
The  locus  of  points  contributing  to  the  power  spectral  density  at  any  specified  frequency  within 
the  doppler  band  describes  a  circle  or  an  arc  of  a  circle,  the  locus  being  the  projection  of  a 
cone  on  that  portion  of  the  range  gate  sphere  lying  within  the  chaff  cloud. 

In  determining  the  doppler  power  spectral  density  function  we  concern  ourselves  with  the 
frequency  band  referenced  to  the  transmitter  frequency  f0.  If  the  maximum  doppler  shift  is 
less  than  f0,  the  return  signal  is  narrow  band  and  two  distinct  bands  of  frequencies,  one  posi¬ 
tive  f.,  the  other  negative  f_,  are  discernible.  Thus,  with  respect  to  the  doppler  frequency 
variable  fjj,  we  may  write  f±  =  ±(f  Q  +  f^). 

The  incremental  power  density  due  to  return  from  an  element  of  arc  lying  along  the  pre¬ 
scribed  locus  of  points  may  be  written 

dP(fD)-^  (fi)G!(0,^)dS,  (65) 

where,  in  spherical  coordinates. 


d  s  =  /sin*  6  d  +  d  92 


The  power  density  from  the  entire  locus  of  equal  doppler  shifted  points  becomes 


I 


P  (fn)  ""  (e)  G2  (5,  <P)  /sin2  9  d  <p2  +  d02. 


(66) 


where  C  defines  the  path  of  integration. 


O 

Wong,  J.  L.,  Reed,  I.  S.,  ond  Koprielion,  Z.  A.,  ’A  Model  for  Rodor  Echoes  from  o  Rondom  Collection  of  Rotating 
Dipole  Scotterers,*  USCEE  Report  159,  Jon.  1966. 
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Defining  #c  as  the  angle  between  position  vectors  to  the  contributing  dipoles  and  the  radar 
velocity  vector,  we  observe 


where 

fD  is  the  doppler  ^  ^aency  shift, 

f0  is  the  transmitted  frequency  (center  frequency  for  a  narrow  band  signal),  and 
V  is  the  magnitude  of  the  radar  velocity. 

The  locus  of  points  C,  defining  the  integration  path,  may  be  found  by  setting  the  normalized 
inner  product  of  the  contributing  dipole  position  vectors  and  the  radar  velocity  vector  equal  to 
cos  0Q.  Hence,  we  obtain 

cos  =  (cos  &Q  +  cos  9  cos  T^/sin  9  sin  y ,  (68) 

cos  y  +  cos  #„  cos  9  „ 

d  <p  = - 5 -  [sin2  0  sin2  y  -  (cos  9  +  cos  9  cos  y)2]  1/2  69.  (69) 

sin# 


Inserting  equation  (69)  in  equation  (65),  we  have 

c-jj  (#)  G2  (5,  <t> )  sin  9  sin  6Q  d  a 


P  (fn> 


r- 

/sin 


i2  9  -  2  cos  #0  cos  y  -  cos2  y  -  cos2  0Q 


where  9  =  cos  1  (-h/d). 

For  y  =  n/2,  horizontal  flight,  this  becomes 


(70) 


o°  (#)  G2  (9,  sin  9  sin  90  69 

P(fD)~  1  ” 


fw  a0 

J  - 


/sin2  9  -  cos2  9„ 


(71) 


For  y  =  0,  vertical  flight,  we  observe  that  7T-0=0O,  and  obtain 


P  (fD)  ~  (77  -  90)  s 


in#p  f 

•'O 


2  7T 


G2  (w  -  #0,  <£)  d 4>, 


(72) 


where  0  <  6’  <>  cos'1  (h/d). 
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For  G  independent  of  <£,  as  in  the  cases  previously  introduced,  we  have 


P  (f0>  i77  ~  60 >  Sin  90  °2  (77  " 


This  becomes,  for  a  uniform  density  chaff  cloud  with  uniformly  distributed  dipole  orienta¬ 
tions,  with  the  range  gate  completely  submerged,  and  with  an  isotropic  antenna. 


P  (fD)  ~  si 


sin  L-‘ 

L  \2Vfo/J 

A-  » 

V  4  V2  f  2 


where  f  is  the  maximum  doppler  shift. 

It  should  be  observed  that  the  calculated  chaff  return  power  spectral  density  is  an  even 
function  with  identical  doppler  spectra  appearing  as  mirror  images  at  both  +f0  and  -fQ  on  the 
frequency  axis.  Considering  the  case  of  vertical  flight,  we  have  plotted  the  doppler  spectrum 
for  an  isotropic  radar  antenna,  Rayleigh  chaff  cloud,  and  uniform  dipole  orientation  distribu¬ 
tion  in  figure  28.  Figure  29  is  a  similar  plot  for  a  cos2  6  antenna  pattern,  Rayleigh  chaff 
cloud,  and  uniform  dipole  orientation  distribution. 


NORMALIZED  DOPPLER  FREQUENCY  l0''0M 

Figure  28.  Doppler  spectral  density  function  for  range-gated  response  at  d  =  100  m  with  an  iso¬ 
tropic  antenna  patt-  n  (Rayleigh  cloud  model,  a  =  0.001)  uniformly  distributeddipoleorienations. 


It  is  relevant  to  the  problem  of  chaff  detection  to  observe  that  the  return  signal  to  the 
radar  will  have  the  maximum  doppler  frequency  shift  and  minimum  frequency  dispersion  when 
the  radar  first  encounters  the  chaff  cloud.  At  this  time,  the  envelope  coefficient  of  correla¬ 
tion  from  pulse  to  pulse  will  be  nearly  unity.  After  the  radar  range  gate  is  completely  sub¬ 
merged  in  the  chaff  cloud,  the  frequency  dispersion  will  reach  a  maximum  and  the  pulse  to 
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NORMALIZED  DOPPLER  FREQUENCY  f0/fDU 

Figure  29.  Doppler  spectrol  density  functior.  for  rc..  e- gated  response  at  d  =  100  m  with  a  cos^d 
antenna  pattern  (Rayleigh  cloud  model,  a  =  0.001)  uniformly  distributed  dipole  orientations 


pulse  envelope  correlation  will  reach  a  minimum.  For  a  radar  with  a  high  loop  sensitivity 
relative  to  the  chaff  return  loop  attenuation,  the  initial  detection  probability  will  be  lower 
than  that  obtainable  from  independent  pulse  returns.  Conversely,  for  a  radar  with  low  loop 
sensitivity  relative  to  the  chaff  return  loop  attenuation,  the  initial  detection  probability  will  be 
higher  than  that  obtainable  from  independent  pulses. 


7.  CORRELATION  FUNCTION 

To  calculate  second-order  probabilities  for  chaff  return  powers,  we  must  first  obtain  the 
coefficient  of  correlation  at  the  sampling  interval.  Since  the  correlation  function  for  predetec¬ 
tion  powers  is  the  same,  except  for  a  scale  factor,  as  the  correlation  function  for  voltages  after 
square-law  detection,  we  obtain  both  results  simultaneously.  Recall  that  the  correlation  func¬ 
tion  of  the  output  of  a  full-wave,  square-law  device  with  a  gaussian  input  becomes10 


Ry  (T)  "•  +  2  Rjj  (t). 


(75) 


where  Rx  (t)  is  the  input  correlation  function  and  a'- j  =  Rx  (0).  The  output  coefficient  of 
correlation  becomes 


Ry  (r>  ~  (®)  _  R*  (7> 

Ry  (0)  -Ry  '  =)  =  TT 


(76) 


^Davenport,  W.  8.  and  riot,  W.  L.,  An  Introduction  to  the  Theory  of  Rondo  n  Signals  end  Noise,  McGraw-Hill,  New 
York,  1958. 
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which  may  be  rewritten  by  using  the  Weiner-Khintchine  relations  as 


Py  CT) 


£ 


(f ')  Px  (f  -  f')  e27,iT(  d  f'  d  f. 


(77) 


That  is,  the  coefficient  of  correlation  for  the  output  of  a  square-law  device  is  the  normal¬ 
ized  Fourier  transform  of  the  self-convolved  power  spectral  density  for  the  bacKscattered 
signal.  The  coefficient  of  correlation  for  the  output  of  a  square-law  device  followed  by  an  ideal 
low-pass  filter  is  the  coefficient  of  correlation  for  sampled  return  powers  and  becomes 


P  (r) 


Px  (f')Px  (f-  f')  e2”iT{  df'  df. 


(78) 


Thus,  from  a  knowledge  of  the  doppler  dispersion  spectral  density  developed  in  the  previous 
section,  the  use  of  equation  (78)  allows  us  to  readily  calculate  the  return-power  coefficient  of 
correlation. 

In  the  case  of  the  uniform  density  chaff  cloud  with  uniformly  distributed  dipole  orientations, 
with  the  range  gate  completely  submerged,  and  with  an  isotropic  antenna  we  may  apply  equation 
(78)  directly  to  obtain 


P(T)  = 


J?(2Wp,7) 
(n  rf 


(79) 


where  J,  is  the  bessel  function  of  the  first  kind  of  order  one,  and  p  (r)  is  the  coefficient  of 
correlation  for  return-power  samples  separated  by  the  time  interval  r . 

For  the  more  general  cases  involving  nonuniform  cloud  densities  and/or  anisotropic  an¬ 
tenna  patterns,  the  direct  application  of  equation  (78)  will  not  usually  be  possible.  In  such 
cases,  the  use  of  fast  Fourier  transform  (FFT)  numerical  computing  technique  as  described 
by  Cooley  et  al.u  may  be  appropriate.  The  computational  procedure  involves  the  self¬ 
convolution  of  the  return-signal  power  spectral  density,  the  separation  of  the  resulting  spec¬ 
trum  into  low-frequency  and  high-frequency  parts,  and  application  of  the  FFT  to  the  low- 
frequency  part. 

Considering  the  examples  used  in  the  previous  section,  we  have  employed  the  FFT  tech¬ 
nique  and  plotted  the  power  correlation  functions  for  the  radar  height  ranging  from  initial 
entry  of  the  range  cell  into  the  chaff  cloud  to  the  total  submersion  of  the  range  cell  in  the 
cloud.  The  correlation  function  for  a  vertical  flight  path,  uniform  dipole  orientation  distribu¬ 
tion,  Rayleigh  chaff  cloud,  and  isotropic  antenna  is  plotted  as  figure  30.  Figure  31  is  a  similar 
plot  for  a  cos2 6  antenna  pattern.  Both  plots  are  mean  adjusted  to  zero  and  the  variance  at 
each  height  normalized  as  in  figures  23  and  24,  respectively. 


IJCooIey,  J.  IV.,  Lewis,  P.  A.,  and  Welch,  P.  D.,  “Application  of  the  Fast  Fourier  Transform  to  Computation  of 
Fourier  Integrals,  Fourier  Series,  and  Convolution  Integrals,'  IEEE  Trans,  Audio  and  Electroacoustics,  vol. 
AU-15,  No.  2,  June  1967. 
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Figure  30.  Correlation  function  for  range-gated  response  at  d  -  100  m  with  an  isotropic  antenna 
pattern  (Rayleigh  cloud  model,  a  =  0.001)  uniformly  distributed  dipole  orientotions 


Figure  31.  Correlation  func>-on  for  range-gated  response  at  d  =  100  m  with  a  cos20  an  ten  no  pctfern 
(Rayleigh  cloui  model,  a  ■-  0.001}  uniformly  distributed  dipole  or*entafions 


From  figure  30,  we  may  observe  that  approximately  5(c/2V)  rf  periods  must  elapse  be¬ 
tween  essentially  uncorrelated  samples  when  the  radar  is  75  m  above  the  chaff  cloud.  After 
the  radar  has  penetrated  into  the  cloud  a  depth  of  100  m,  however,  approximately  four  essen¬ 
tially  uncorrelated  samples  may  be  obtained  in  the  same  time  interval.  Hence,  if  the  detection 
process  involves  the  incoherent  integration  of  a  number  of  successive  pulses  and  the  pulse 
repetition  period  is  shorter  than  5(c/2V)  rf  periods,  the  fluctuations  about  the  mean  of  the 
integrated  process  will  be  considerably  greater  at  a  height  of  75  m  than  at  -100  m. 


8.  SUMMARY 

In  this  report,  we  have  provided  a  model  for  chaff  sufficiently  general  to  characterize  a 
diverse  set  of  distribution  geometries.  Simultaneously,  we  have  retained  the  requisite  sim¬ 
plicity  of  description  so  that  the  methodology  of  specific  calculations  is  not  unnecessarily  en¬ 
cumbered.  We  have  introduced  and  developed  the  concept  of  the  chaff  impulse  response,  noting 
that  such  a  function  is  theoretically  fictitious,  but  recognizing  it  as  a  practical  approximation 
for  chaff  backs catter  to  a  short-pulse  radar.  A  significant  difference  between  the  delay  time 
to  peak  for  the  chaff  impulse  response  as  compared  with  the  rough  surface  scatter  impulse 
response  has  been  noted. 

We  have  also  examined  the  response  to  a  chaff  cloud  as  viewed  from  a  fixed  range  cell  as 
the  radar  moves  relative  to  the  cloud.  The  probability  density  functions  of  such  a  signal  have 
been  reviewed,  and  the  spectral  density  function  for  this  3ignal  in  consonance  with  our  chaff 
cloud  model  has  been  derived.  Considering  square-law  detection  of  this  signal,  we  have  ob¬ 
tained  the  correlation  function  of  the  detector  output  and  examined  its  fluctuation  rate  as  the 
position  of  the  radar  above  and  within  the  chaff  cloud  is  varied. 
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GLOSSARY  OF  TERMS 

a  =  radius  of  the  circular  cross  section  of  a  thin  wire  chaff  dipole. 

C^  =  proportionality  factor  for  maximum  value  of  the  ratio  5^  A2  • 
c  =  speed  of  light  in  vacuo. 

d  =  range  gate  nominal  distance,  distance  from  the  radar  to  chaff  dipoles  con¬ 
tributing  to  the  return  at  a  specified  time. 

Ad  =  incremental  distance,  depth  of  the  range  resolution  cell. 

E  {•)  =  statib-ioal  expection  operator. 

FFT  =  fast  Fourier  transform  numerical  computing  technique, 
f  =  frequency. 
fD  =  doppler  frequency  shift, 
f  =  maximum  doppler  frequency  shift. 

f  =  frequency  of  transmitted  wave  (center  frequency  for  a  narrow  band  signal). 
G  “  antenna  power  gain. 

G0  -  peak  value  of  antenna  power  gain. 

Ut  receiving  antenna  power  gain. 

Ut  -  transmitting  antenna  power  gain. 

h  =  distance  from  the  radar  to  the  chaff  cloud  top  (positive  when  the  radar  is 
above  the  cloud). 
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h'  =  a  vertical  distant  measured  from  the  radar. 

I  -  length  of  chaff  dipole. 

m  (• )  =  mean  value  function  of  a  stochastic  process. 

n  =  total  number  of  chaff  dipoles  contained  in  the  vertical  column  over  each  hori¬ 
zontal  unit  area. 


N .  =  a  Poisson  distributed  random  variable,  the  number  of  chaff  dipoles  contained 
in  region  v. 

N  (v  )  =  a  Poisson  process  specifying  number  of  chaff  dipoles  contained  in  a  region 
centered  on  point  P. 

Pr  =  backscattered  received  power,  square  of  envelope  of  the  return  signal. 

P*  =  chaff  impulse  response  function. 

P  =  transmitted  power,  the  square  of  envelope  of  the  transmitted  signal. 

Px  ( f )  =  power  spectral  density  function  of  the  random  process,  X(t). 

p  (x)  =  first-order  probability  density  function  of  the  random  variable,  X  [stationary 
random  process  X(t)] 

P2  (Xr  Xj)  =  second-order  probability  density  function  of  the  random  variables  Xj  and  X2 
[random  process  X(tj )  and  X(t2)] . 

R  (X,  V)  =  a  function  of  the  wavelength  \  and  the  polarization  state  r. 


R  (X.  r;)  = 


c  G  l  X4  C, 


4  (4  77)2 

Rx  (r)  =  correlation  function  for  the  stationary  random  process  X(t). 
t  =  time  variable. 

V  *>  magnitude  of  radar  velocity. 


v  =  mean  chaff  'll  speed. 

a  =  a  configuration  parameter  related  to  the  chaff  cloud  density  gradient. 
y  -  flight  angle  between  the  vertical  and  the  radar  velocity  vector. 
h  =  chaff  dipole  volume  density. 

S  m  =  maximum  chaff  dipole  volume  density. 

6  =  elevation  angle  of  a  radar-centered  spherical  reference  frame. 
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&'  =  elevation  angle  of  a  dipole-centered  spherical  reference  frame. 

\  =  nominal  wavelength  (wavelength  associated  with  the  center  frequency  of  a 
narrow  band  signal). 

p  (r)  =  coefficient  of  correlation  of  a  stationary  random  process, 
a  =  monostatic  radar  cross  section. 
g  =  mean  value  of  the  monostatic  radar  cross  section. 

g  =  mean  value  of  the  monostatic  radar  cross  section  for  a  horizontally  polarized 
incident  wave. 

?v  =  mean  value  of  the  monostatic  radar  cross  section  for  a  vertically  polarized 
incident  wave. 

=  the  chaff  volume  radar  cross  section  for  a  v  polarized  incident  wave. 

&x  =  standard  deviation  of  random  variable  X  [stationary  random  process  X(t), 

^=«x<W 

r  =  time  interval  between  successive  samples,  the  pulse  repetition  period,  the 
elapsed  time  since  chaff  dispersal. 

4>  =  azimuth  angle  of  a  radar-centered  spherical  reference  frame. 

4>‘  -  azimuth  angle  of  a  dipole-centered  spherical  reference  frame. 

tp  =  a  uniformly  distributed  (0,  2tj)  random  variable,  a  random  angle  of  rotation 
measured  in  a  horizontal  plane. 
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APPENDIX  A.  IMPULSE  RESPONSE  FUNCTIONS 

Following  are  integrated  egressions  for  the  impulse  response  functions  developed  in 
section  4  and  plotted  as  figures  10  through  17  (body  of  report). 

A-l.  Exponential  Cloud  Model  S  =  &a  (1  -  e-a<h  ~h)) 

Following  are  mathematical  derivations  for  uniform  distribution  of  dipole  orientations  or 
horizontal  distribution  of  dipole  orientations  with  horizontally  polarized  wave 

A-1.1  Isotropic  Antenna  Gj  =  G10 


A-1.5  Horizontal  Distribution  of  Dipole  Orientations  with  Vertically  Polarized  Wave  — 
Isotropic  Antenna 


“i  =  Oi, 
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A-2  Rayleigh  Cloud  Model  8  =  na  (h'  -  h)e"°/2  <h'-h> 

Following  are  mathematical  deviations  for  uniform  distribution  of  dipole  orientations  or 
horizontal  distribution  of  dipole  orientations  with  horizontally  polarized  wave 

A-2.1  Isotropic  Antenna  Gi  =  c,io 

P*  (H)  =  R  (K-  lH  (1  -  e-a/2<d-h)*]f  -dih^h 

d3 

.  R  V)  [e-a/2  (dth)2  +  e-o/2  (d-h)2].  h  <  -  d 
d3 

=  0.  h  >  d 


A-2.2  Double- Lobe  Antenna  G,  =  G20  cos2  6. 


>s  (d)  g_Lx:  VI  jjh4  +  12J?°  +  -  jd4  *  1  (d2  *  h  d  +  h2> 


-  d  S  h  S  d 
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A-3.  Conclusion 


Although  the  exponential  chaff  cloud  model  leads  to  an  inverse  square  range  attenuation 
law  in  the  limit  h--a>,  the  Rayleigh  chaff  cloud  model  is  bounded  above  by  adverse  cube 

zss. “•  re8po*1  *• 


